#!/usr/bin/env python 

import numpy

from specks import *

s = Speck()

for i in range(13):
    s.addAtom('Al', numpy.random.normal(0, 2, 3))

xyz.save("qm.xyz", s, 'w')

beam.speck = s
beam.start()

m = Morse()

qm = QuickMin(s, m)

m.force(s)
while norm(s.f) > 0.0001:
    qm.step()
    print norm(s.f)
    xyz.save("qm.xyz", s, 'a')
    



